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Abstract 

Several aspects of interplay between Anderson localization and coherent am- 
plification/absorption, and aspect of mirrorless laser for a laser-active (ampli- 
fying) disordered dielectric medium have been addressed. We have calculated 
the statistics of the reflection coefficient and the associated phase for a light 
wave reflected from a coherently amplifying/absorbing one-dimensional disor- 
dered dielectric medium for different lengths, different disorder strengths and 
different " amplification" /" absorption" parameters of the disordered sample. 
Problems with modeling coherent amplification/absorption by complex po- 
tential have been discussed. We have shown an alternative way to model the 
coherent absorption. Several conceptual and physical aspect for coherently 
amplifying/absorbing media with disorder also have been discussed. 



Typeset using REVTgX 



I. INTRODUCTION 



Wave propagation in disordered media and its consequences to the disorder induced lo- 
calization, or Anderson localization have been studied in detail for electron (quantum) 
waves 0, as well as for light (classical) waves for last four decades. We have discussed 
several aspects of transport and localization for ID electronic systems in our previous paper 
In this paper we are dealing with an interesting aspect of the application of the An- 
derson localization to enhance the coherent light amplification in a coherently amplifying 
disordered medium. Due to the bosonic nature of the light quanta, there is possibility of 
coherent amplification or coherent absorption(attenuation) for light; and these possibilities 
are absent for electronic (fermionic) case. Background to study this problem comes from 
the following two facts: (1) coherent amplification in an coherently amplifying medium is 
a non conserving scattering process where temporal phase coherence of a wave is preserved 
despite amplification and (2) coherent back scattering (CBS) in a disordered media is the 
main cause of weak and strong localization will not be affected with additional presence of 
a coherently amplifying medium due to the persistence of phase coherence of the interfering 
waves. The preservation of the CBS, hence the localization, despite amplification can bring 
the interesting possibility of obtaining synergitic enhancement of wave amplification , or 
laser action without mirrors due to the confinement by the Anderson localization, in an 
optically pumped laser active disordered media. In our previous short paper p such possi- 
bility was shown theoretically for ID amplifying disordered media. Our calculation lead 
to an important result that in the presence of amplification in a disordered optical medium, 
the Anderson localization enhances the coherent amplification. In fact, recent experiments 
||6HTT| support aspects of the lasing action, or more precisely, enhanced Amplification of 
Spontaneous Emission (ASE), in an amplifying and strongly scattering medium without an 
external resonant cavity. There is current interest to see the affect of disorder on coher- 
ent amplification/absorption in different disorder regimes. In this detailed paper we have 
discussed the possible background and probable several physical aspects of the problem of 
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coherent light amphfication using the Anderson locahzation and a detailed study of the ID 
amplifying/ absorbing disordered medium. 

It is now well known that in ID and 2D all states of a wave are localized for both elec- 
tronic and optical disordered systems 0,0]. This motivate one to study synergetic effect of 
amplification and localization for lower dimensions. We study analytically and numerically, 
the probability distribution of the reflection coefficient (|-Rp) and phase (6) of the complex 
amplitude reflection coefficient {R = |i?|e*^) of a wave reflected from a one-dimensional 
coherently absorbing/attenuating disordered medium for a Gaussian white-noise potential 
for different lengths of the sample, with different strengths of disorder, and with different 
active (gain or loss) parameter values. We set up a general framework to study ID coher- 
ently amplifying/ absorbing disordered media and discuss the context of physical realization 
for both optical and electronic cases. Coherent ampliflcation we mean mainly is due to 
the stimulated emission of radiation which can only be applicable to an optical (bosonic) 
case and has no electronic counterpart. Coherent absorption can have physical situation 
in both the optical and the electronic cases — for the electronic(fermionic) case it signifles 
effectively the incoherent part of the electron wave function caused by inelastic scattering. 
We consider here mainly the ID Helmholtz equation, where the coherent (linear) ampli- 
flcation/absorption is modeled by adding a constant imaginary part to the real potential. 
We derive flrst a Langevin equation for the complex amplitude reflection coefficient R{L) 
and then the Fokker-Planck equation in the (|-Rp, 6')-space for evolution with the length 
of the sample for different parameter values: (1) disorder strength, and (2) active parame- 
ter. We discuss the appropriate physical situations for both the cases, i.e., absorption and 
amplification. 

Modeling amplification/absorption by a complex potential, however, always gives a re- 
flection part due to the mismatch between the real and the imaginary potentials. We discuss 
these issues in detail. We also derive a Langevin equation and then the FP equation by a 
phenomenological modeling where absorption has no concomitant reflection part. We also 
have discussed recent experimental and theoretical works of others on coherent ampliflca- 
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tion / absorption. 

II. COHERENT BACKSCATTERING IN THE PRESENCE OF COHERENT 
AMPLIFICATION/ABSORPTION IN AN ACTIVE DISORDERED MEDIUM 

The main cause of the weak and the strong locahzation in a disordered medium is 
coherent-back-scattering (CBS) [[T2i p|. In the weakly locahzed regime, CBS is the correction 
to the classical transport, and in the strongly localized regime CBS is the dominant effect. 
In the presence of coherent amplification/attenuation, the partial wave amplitudes get am- 
plified/attenuated in the same way along the echo (i.e., time-reversed ) paths when they 
meet at the starting point tq as shown in Fig.|l|, i.e., the phase is not lost due to the coher- 
ent nature of the amplification, or attenuation which is relevant to the stimulated emission 
of radiation, or the coherent (stochastic) absorption. Let now briefly describe disordered 
coherently amplifying/absorbing media relevant to this problem. 

HI. AMPLIFICATION 

A. Light localization and coherently amplifying media 

By coherently amplifying media we mean stimulated emission of radiation in a lasing 
medium. Here, we show that the presence of localization in a highly scattering medium, 
when the medium is also an active amplifying one, can enhance the amplification. There is 
the possibility of using Anderson localization as a " non- resonant feedback" mechanism for 
lasing action. Recent experiments which show self-sustaining lasing action in an optically 
amplifying and highly scattering medium, without an external cavity, support these aspects 
which we will discuss in later section. Although the real situation of stimulated radiation in 
a lasing medium involves non-linear amplification, for the sake of simplicity we will consider 
a linear amplifying medium only and discuss enhanced amplification by the disorder induced 
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confinement. By the linear medium here we mean the refractive index being independent of 
the light intensity. 

Light localization is an interesting and clean problem due to the fact that the energy 
of the photon is high and no effect of temperature and photons are noninteracting. There 
are, however, problems in localizing light 0|T3. This is because the "effective potential" 
of a dielectric medium depends on the incoming wave energy, or the wavelength (A). The 
problem is to satisfy the Mott-Ioffe-Regel (MIR) criterion for the localization which is 
essential for any wave localization. The MIR criterion for the localization demands that 
kle << 1, where k is the wave number of the incident wave and le is the elastic (coherent) 
scattering mean free path for transport as shown in Fig.^ . For longer wavelengths, or in the 
low energy limit of the incident wave, scattering is dominated by the Rayleigh scattering, 
Ig ~ A"^"*"^, for a d dimensional scatterer. In the other extreme limit of high frequency, 
or short wave length, the system is almost transparent (non-scattering) to the incoming 
wave for wave lengths smaller than the correlation length of the random potential when 
the response of the system is that of geometrical optics. Localization is, thus, missed in 
both the wave length limits, the large as well as the small. While in ID and 2D all states 
are localized for arbitrary weak disorder, there is a very wide localization window where the 
localization length is very large. In 3D there may be a small free-photon localization window 
for intermediate wavelengths. Here the MIR criterion can be satisfied in an indirect way by 
first creating a pseudo-band gap by fulfilling the Bragg condition k.G = ^G.G through a 
periodic lattice of dielectric medium and then deliberately introducing a weak disorder to 
the system [|15|,|16|. It was shown for 3D that in such systems, the density of states near the 
band gap correspond to localized states. 



B. The questions we are interested 

What will happen to coherent light amplification in an amplifying (active) 
medium, when the light modes are already disorder-localized in that medium? 
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For ID and 2D all modes are localized for light wave and one can see the direct affect of 
localization on coherent amplification/absorption. In 3D there is a small localization window 
and possibility of a physical situation can be as follows. Suppose the wavelength Al of a 
lasing light lies within the localization window. A typical localization window for 3D is 
shown in Fig.|^. Now, consider a three-level lasing scheme as shown in Fig.|^. For lasing 
action electrons have to be pumped first from the ground level to the upper excited levels. 
The transitions from the excited level to the intermediate meta-stable level are arranged to 
be non-radiative. Lasing action occurs between the meta-stable and the ground states. Now 
it is possible to arrange such that the shorter pumping wavelength Xp that excites electrons 
from the ground state to the upper excited states lies within the delocalized geometric- 
optical regime, while the meta-stable to the ground state radiation wave length, i.e. the 
lasing wavelength A^, lies within the localization window. Then, the question is how the 
localization and the coherent amplification synergetically lead to enhanced amplification. 
Also, how does the lasing wavelength get selected, if at all in the absence of a resonant 
cavity structure providing the feedback. 



C. Different feedback mechanisms for lasing action 

In a typical amplifying medium, a wave has to traverse a long path-length for amplifi- 
cation, by stimulated emission. Simultaneously, the wave may lose its energy by other loss 
mechanisms like non-radiative dissipation, escape out of the system etc.. The wave requires 
a positive feedback to gain a higher amplification. Such a feedback can be given in several 
ways. 



1. Resonant feedback 

In the case of resonant feedback, the laser active material is kept inside a resonant cavity 
r?| , which is generally made of two parallel mirrors as shown in Fig.^. Stimulated emission 



has a typical width for its intensity distribution (gain profile) over the wavelengths. Resonant 



cavity is made by matching the spacing between the two mirrors such that it can give positive 
feedback selectively to a set wavelengths that "fit" into the cavity. Given a cavity of length 
L, we will get the resonant positive feedback, if mA/2 = L is satisfied, where m is an integer. 
The wavelength which satisfies the above resonance condition will survive, because only this 
wavelength will overcome the energy losses in the medium through multiple passes. Intensity 
for the rest of the wavelengths will die down as there is no feedback to them to compensate 
for the loss in the medium. The line-width of the laser light will now be determined by the 
cavity — Q(quality)-factor. It is much narrower than the width of the gain profile of the 
laser material. It is also narrower than the widths of the cavity resonances. It is in fact 
" gain-narrowed" . 



2. Non-resonant feedback/feed- forward 

A distributed "gain" to a wave propagating in an amplifying medium can be given to 
it in its forward journey (single pass) also, if it moves over a sufficiently long path length 
as shown in Fig.^. An example of this type of medium is to be found commonly in our 
galaxy. Astrophysical "maser" [l^, like the SiO maser in our galaxy gets feedback by this 



mechanism. The gain produces its own spectral narrowing without any resonant structures. 
In order to see how this happens, note that the medium amplification has a gain profile with 
respect to the wavelength. A starting intensity profile, which also has a distribution over 
the wavelengths, passing trough the amplifying medium is modulated by the gain profile. 
In the process of a long-length journey, only the prominent peak of the initial intensity 
profile will be effectively amplified. Thus, let /(A) be the starting intensity distribution 
of a stimulated radiation and ^'(A) be the gain profile of the medium. After n successive 
coherently amplifying encounters, the intensity distribution will be amplified iteratively as: 
/(A) ^ I{\)g{\)^I{X)g{\r I{\)g{\r. 



The latter would tend to a single delta-function as n — cxd, even when the initial profile 
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/(A) is broad, but singly peaked. 

Fig.^b shows a schematic typical "gain narrowing" by such a "feed- forward" mechanism. 
"Feed-forward" is a better nomenclature here than "feedback". 

A variation on the above feed-forward theme is readily realized in a diffusing, amplifying 
medium where random multiple scattering (turbidity) can enhance the path length traversed 
by light ( or increase its residence time) before emerging out of the amplifying medium. 
Thus if la is gain length for the pure amplifying medium (i.e., the length over which the 
intensity is amplified by a factor e"*"^), and le is the transport mean free path for diffusion 
due to scattering, then for the three-dimensional sample size L > Lsr ~ n^Jlgle/S, the gain 
exceeds the loss at the boundaries and self-sustaining emission — Super Radiance(SR) - can 
result. There is really no feed-back loop here and hence no true laser action. What one has 
is the Amplification of uncorrelated Spontaneous Emission (ASE). It does, however, show 
gain narrowing as discussed above. 

3. Non-resonant feedback by Anderson localization : Mirrorless Random Laser 

We are trying to explore the possibility of using the Anderson localized states for a feed- 
back mechanism in the disordered lasing media. When a wave passes through a disordered 
amplifying medium it undergoes coherent multiple scattering. It is reasonable to think that 
due to the multiple scattering processes a wave can get amplified before it escapes from 
the system. In a strong disordered media, coherent back scattering ensures the wave to 
bring back to its starting point. Now the questions are : (a) is this amplification sufficient 
to overcome the loss in the system, (b) can we use the system (amplifying medium with 
disorder) as a pure optical amplifier, and (c) is there a mechanism for wavelength selection. 
The answers to (a) and (b) seem to be yes, as supported by recent theories and experiments 
and answer to (c) is still unknown. Recent series of experiments by Lawandy et. al. P,|8|JTT[| 
have shown evidence for self-sustaining lasing action in an amplifying, optical scattering 
medium. Though it seems more like Amplified Spontaneous Emission (ASE) |^|TU[. The 
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experiment was performed on rhodeamine 640 perchlorate laser dye, dissolved in methanol, 
as an active medium. Colloidal suspension of Ti02 nanoparticles (coated with AI2O3) in the 
active methanol solution was used as the optical scattering random medium. The emission 
from such a system was found to exhibit multi-mode laser oscillations in the absence of any 
external resonant cavity. With no Ti02 (no strong disorder) in the active solution, there 
was no lasing action from the pure laser active dye dissolved in methanol. 

Very recent experiment by Wiersma et. al. ^ shows enhancement of amplification due 
to coherent back scattering from an coherently amplifying optical disordered solid state 
medium. This experiment is more relevant to our calculation, which deals with amplification 
in a medium due to localization. 

These experiments motivate one to study the phenomenon of amplification-enhancement 
by disorder in more detail. The model we are going to consider is for the case of a pure 
ID random amplifying medium. A physical realization may be a single longitudinal-mode 
optical fiber doped with a laser active material {Er^~^, say) and disordered intentionally. 



IV. COHERENTLY ABSORBING MEDIA 

There have been recent interest in the role of absorption on localization, and on the 
different length scales in the problem in the presence of absorption. For the bosonic case, like 
light etc., a coherent state (e.g. a laser beam) is an eigenstate of the annihilation operator. 
Removal of a photon(absorption) does not destroy the phase coherence. Fermions cannot 
be annihilated in the context we are considering here. For the fermionic case, the physical 
picturization of the absorption will be some kind of inelastic process (like scattering by 
phonons), where electrons lose its partial temporal phase memory. This type of absorption 
is called stochastic absorption. Absorption can be of an other type also, the so called 
deterministic absorption. In this case, particles are shut out by a chopper from its beam 
path. It has been shown that the stochastic absorption retains more coherence than the 
deterministic absorption while absorbing (removing) the same amount of particles . 
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Recent theoretical studies [pO|-p2| have shown that the absorption in case of hght wave 



does not give any cut-off length scale for the localization problem. If a system is in the 
localized state, absorption will not kill the localization to make the system again diffusive. 
A sharp mobility edge exist even in the presence of significant absorption in 3D. 

V. MODEL 

A. Modeling coherent amplification/absorption by complex potential 

Adding a constant imaginary part with the proper sign to the real potential of the 
Maxwell/Schrodinger wave can model the linear coherent amplification/absorption. The 
Maxwell equation can be transformed to the Schrodinger equation where the effective po- 
tential depends on the frequency of the incoming wave. Also, both the Maxwell and the 
Schrodinger equations can be transformed to the Helmholtz equation. 

Adding a constant imaginary potential iVa to the original real potential V{x), the 
Schrodinger equation becomes: 

where fi^ /2m = 1. 

Similarly adding an imaginary constant term ze^ to the original random refractive index 
eo + e{x), the Maxwell equation becomes: 

^ + ^h + eix)+^e^]E =0 . (2) 

Here the subscript 'a' denotes amplification/absorption(attenuation). 
The corresponding Helmholtz equation can be written as: 



+ k'[l + {r]{x) + ir]a)]u = , (3) 



where ri{x)=—V{x)/k'^, ria=—Va/k'^ for the Schrodinger equation, and ri{x)= Va= 
and k"^ = for the Maxwell equation. Here, eg is the constant dielectric background and 
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e(x) is the randomly spatially fluctuating part of the dielectric constant, k = 27r/A, where 
A = wavelength in the average medium (eo). 

The model we are considering here is for one-channel scalar wave where the polarization 
aspects have been ignored. This holds good for a single-mode polarization maintaining 
optical flber. Formally, we will consider ID one-channel Helmholtz equation only, and will 
specify at the proper place which situation is appropriate, optical or electronic. 

It has been pointed out by Rubio and Kumar that modeling absorption by complex 
potential will always give a concomitant reflection. Increasing the strength of the complex 
potential will not increase the absorption. There is a competition between the ampliflca- 
tion/absorption and the reflection. For very high absorption, the system may try to act as 
a perfect reflector. Later we present a detailed analysis for the modeling of absorption by 
"fake", or "side", channels obviating the need for a complex potential. 

B. The Langevin Equation for the complex amphtude reflection coeflicient R(L) 

The Langevin equation for the complex amplitude reflection coefficient R{L) can be 
derived by the invariant imbedding method from the Helmholtz equation Eq.§, 



= 2^kR{L) + ^^(l + r^(L) + ^r^,){l + R{L)f , (4) 
with the initial condition R(L)=0 for L=0. 

C. The Fokker-Planck Equation 

Now, taking R{L) = a/tc*^, the Langevin equation Eq.^ reduces to two coupled differen- 
tial equations. 

dv 1 

— = kr]{L)r^{l — r) sinO — 2kriar 



krjaV^ [1 -\- r) cos 6 (5) 
k 

2k + -1 
dL 2 



^ = 2k + K{L) \2 + cose{r^l^ + r""^ 
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(6) 



The same way as described in Ref. [0], using stochastic Liouville equation for the evolution 
of probabihty density and then integrating out the the stochastic part by Novikov's theorem, 
we get the Fokker-Planck equation: 



9P(r, 



dl 



sin^|-(rV2(l-r)) + | 



-\ 2 



Pfr, 



D 



d 



(7) 



where we have introduced the dimensionless length / = j and D 



qk 



i/ia, and i 



i^qk^) ^ is the localization length and = -^t^ is the active (amplifying/absorbing) length 
scale in the problem. 



D. Parameter scales in the problem 

The above Fokker-Planck equation (Eq|^ has effectively three dimensionless scales: 

(1) Length I = L/C^, where L is the length of the sample and ^ = {■^qk'^)~^ is the localization 
length, k is the incoming wave vector and q is the delta-correlation strength of the potential. 

(2) Disorder parameter 2k^. 

(3) Active parameter D = ^/C,a, where C,a = (2a^ active length scale. 

In the case of amplification, C,a is usually known as the gain length Ig. The gain length Ig 
decreases with increasing population inversion (pumping), and eventually saturates at high 
enough pumping rate. 
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VI. RESULTS AND DISCUSSION 



A. Analytical solution of the Fokker-Planck equation within the random phase 

approximation (RPA) 

The Fokker-Planck equation Eq.^ can be solved analytically for the asymptotic limit of 
large lengths in the random phase approximation (i.e., for the weak disorder case). Inte- 
grating the phase part of the FP Eq.|^, i.e., P(r) = 1/ {2n) J P{r, 9)d9, one gets FP equation 
of the reflection coefficient r: 

dl ^ ' dr^ 

+ [l + (-6 + D)r + 5r2]^^^ 

or 

+ [{-2 + D)+4r]P{r,l) , (8) 

(To make the point more clear we have written P(r) as P(r, /) to show the explicit length- 
dependence) 

Here, D < corresponds to coherent amplification, D > corresponds to coherent 
absorption, and D = corresponds to the unitary case. To emphasize the dependence of 
P(r, I) on D we will write here P(r, /) = P^{r, I) . 

Few asymptotic features of P^{r,l) can be obtained directly from Eq.^). We can note 
that Eq. (H) reduces to the unitary case, for {D = 0) in the limits rja = trivially, and 
large q nontrivially. The latter case is dominated by strong disorder and the localization 
length is too short for the wave to have penetrated the amplifying medium. Then, the 
statistics P^{r,l) have steady state distribution for / —>■ oo. It can be obtained by setting 
dP{r, l)/dl = and solving the resulting equation analytically. We get the following steady 
state distributions: 
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[a) Amplifying media : 



lim P^(r,/) 

l—*oo 
D<0 



\D\exp{' 



(l-r)2 





for r > 1 
for r < 1, 



(9a) 



(6) Absorbing media : 



lim P^(r,/) 



/ — fOO 

D>0 



I — >oo 
D=0 



\D\exp{\D\)exp{--^) 





for r < 1 
for r > 1, 



l^c) Passive media : 

lim P^(r,/) =5(l-r) 



(9b) 



(9c) 



In Fig.p we have plotted P^{r,l = oo) for the three limiting cases with parameter 
D = +1, 0, —1. These steady state probability distributions for P{r) with the weak disorder 
show that for the case of amplifying media, the probability distribution of the reflection 
coefficient effectively lies within r > 1 region and spread out to infinity. In this asymptotic 
limit of large sample length the average < r > diverges for any D, i.e, we have superradiance 
in fact laser action. 

Fig.^ clearly shows that the Anderson localization enhances the amplification. Thus, 
there is a possibility that the Anderson localization can be used for the non-resonant feedback 
mechanism for lasing action. These probability distributions have universal nature and have 



been studied recently in more detail using different approach by several authors p5|-[3q . 
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B. Numerical solution of the Fokker-Planck equation beyond the 
random-phase-approximation: Marginal probability distribution P{9) of the phase 6 
and the probability distribution P{r) of the reflection coefficient r. 

The Fokker-Planck equation, Eq.|^, is difficult to solve analytically beyond the random 
phase approximation. We solve the full Fokker-Planck equation (|^ numerically to inves- 
tigate the effects of phase-correlation on the statistics of the reflection coefficient. Initial 
probability distribution P{r, 9) at I = for our numerical work was taken to be the same 
as in the case of the passive disordered media described in our previous paper [Q. Now, 
r in case of absorption is limited between and 1, and in case of amplification, r is lim- 
ited between and oo. For our numerical calculation for the amplifying case, however, we 
have considered values of r sufficiently large but finite such that the truncated range of the 
probability distribution shows the main feature of the full distribution. Below we are giving 
in detail, the simulation results of the Fokker-Planck equation in all important parameter 
regimes. 

1. Phase distribution with active parameter D, for the case of weak disorder 

In Figs.|^(a),(b) we have plotted the variation of the phase in the random phase ap- 
proximation with the active parameter with fixed disorder parameter 2kE, = 50 (weak 
disorder regime) and sample length /. It is clear that the random phase approximation holds 
up to a large value of \D\, for both (a) amplifying and (b) absorbing medium. At this point 
we can tell that the demand of the Ref. [^,^] on the issue of phase distribution is invalid. 
Once disorder is weak, the RPA persist to a higher amplification/absorption values. 

2. Weak disorder with weak amplification/absorption 

Figs.^ {a.),{h) are the plots of the reflection coefficient r distribution P(r) for weak 
l^l. These plots show the evolution with the length of the sample and saturation for large 
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lengths. The distribution matches the analytical result for the asymptotic limit of large 
sample lengths within the random phase approximation for both (a) the amplifying and (b) 
the absorbing cases. 

3. Weak disorder with strong amplification/absorption 

Figs.P(a),(b), are the plots of P(r) for different lengths, (a) In the case of an amplifying 
medium P(r) has complicated behavior, (b) In the case of absorbing media, the waves get 
absorbed before it gets reflected which explain the peak of the P(r) near r = 0. 

4- Strong disorder with weak amplification/absorption 

In Figs.p!0|(a),(b), two top plots are the evolution of P(r) with different sample lengths 
/ for (a) Amplification and (b) Absorption. Competition between the localization and am- 
plification/absorption is implied by the double peaked distribution. The signature of more 
localization is that the reflection probability will peak near r = 1 and the signature of more 
absorption is that P(r) will peak near r = 0. 

In Figs.|TO|(a),(b), bottom two plots are the phase distributions P{9) correspond to the 
P(r) distribution for (a) amplifying and (b) absorbing cases. The P{0) distributions are 
essentially like those of the passive medium ^ , peak near 9 = n symmetrically and do not 
get distorted much. Cases are for: (a) Amplification and (b) Attenuation. 

5. Strong disorder with strong amplification/absorption 

In Figs.^l](a),(b), top two plots have the following features: 

(a) for the amplifying case the probability spreads to the right side of r = 1 for / = 1 and 
peak of the probability P(r) slowly moves towards r = 1 for larger lengths. 

(b) shows that for the absorbing media, the probability initially peaks near r = and the 
peak slowly moves towards r = 1 for larger lengths. 
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In Figs.|TT|(a),(b), bottom two plots show that the steady state phase distributions is a 
delta function peaking symmetrically about 6 = tt, for both the cases (a) amplifying and 



(b) absorbing. This result contradict the recent results of Ref. [p9| ,p0|] 

These results (Figs. 0(a),(b)) support that the system is behaving as a perfect reflector. 
It was indicated by Rubio and Kumar PB[ , that modeling with a complex potential is always 
associated with reflection due to the mismatch of potentials (real and complex). Increasing 
the strength of the imaginary part of the potential will not necessarily increase the amount 
of absorption or amplification, i.e., absorption/amplification is not a monotonic function of 
the strength of the complex potential. Here we show by explicit numerical simulations that 
the system behaving as a perfect reflector for higher lengths, with strong disorder and strong 
active parameter \D\. We will discuss this issue next. 



VII. STOCHASTIC ABSORPTION 

A. Stochastic Absorption from an Absorbing Side ("Fake") Channel (modeUng 

absorption without reflection) 

We have shown and discussed that the absorption is not possible without reflection in a 
model calculation where amplification/absorption is modeled by adding a constant imaginary 
potential to the real potential. In the Langevin equation for R{L), derived from these types 
of model, the wave always gets a reflected part along with the absorption. There is, however, 
another way of deriving a Langevin equation for the reflection amplitude R{L) such that the 
absorption does not have a concomitant reflected part. This approach is motivated by the 



work of Buttiker p9| , ^ where some purely absorptive "fake" (side) channels are added to 
the purely elastic scattering channels of interest. A particle that once enters to the absorbing 
channel, never comes back and it is physically lost. In the Appendix, we have derived the 
Langevin equation for R{L) following the Buttiker approach. This Langevin equation has 
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some formal differences from the Eq.^. 



B. The Langevin Equation 



The Langevin equation for the stochastic absorption only turns out to be : 

= -aR{L) + 2ikR{L) + i^HL)){l + R{L)f , (10) 

(for derivation and discussion see Appendix) with the initial condition R(L)=0 for L=0, and 
a is the absorbing parameter. 



C. The Fokker-Planck Equation 

Similarly, following the same steps as we discussed previously, from the Langevin Eq(p!OD 
we derive the Fokker-Planck equation: 



dP{r, I 



P(r, 



where parameters / and 2k^ are same as defined in Eq.|^, and D = 



D. Solution of the FP (Stochastic absorption case) equation within the Random 

Phase Approximation 

Eq |ll] gives the same Fokker-Planck equation in r, i.e., Eq. ^ for the reflection coefficient 
when phase part is integrated out by the RPA in the limit of weak disorder. Hence, steady 
state solutions are also same as Eq.pl]. 



18 



E. Numerical solution of the FP equation for the case of stochastic absorption 

beyond the RPA for strong disorder 

For the weak disorder case, where random phase approximation is vahd, results are 
the same as for modehng with a complex potential. We will consider here only the strong 
disorder case for the numerical calculations of the Fokker-Planck equation ^ within the RPA. 

1. Strong disorder and Weak Absorption 

Fig.^(a) is the plots of P(r) and P{9) , respectively, for the strongly disordered media 
with weak stochastic absorption. Probability distribution for the reflection coefficient has 
double peaked behaveour which arises due to the competition between the absorption and 
localization. 

Fig.|^(b) is the plot of P{9) which is typical of a double peaked symmetric distribution, 
similar as that for strongly absorbing and strong disordered media. 



2. Strong disorder and Strong Absorption 

Fig.^(a) is the plot of P(r). The wave is absorbed in the medium before it gets reflected 
as we have modeled absorption without reflections and we are considering the case of strong 
absorption, and hence the probability of absorption is more than that of reflection. 

Fig.0(b) shows the phase distribution P{6), which is a double peaked symmetric distri- 
bution and differs from the model of absorption by imaginary potential, which was a delta 
function distribution ai 6 = n. 



VIII. DISCUSSION AND CONCLUSIONS 

There are a few recent theoretical works on random amplifying media for the weak and 
the strong disordered regimes. Zhang p5|-p^ has simulated the ID disordered amplify- 
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ing/absorbing medium by transfer matrix method and got the same steady state distribu- 
tions for P(r) as ours, thus confirming our analytical treatment. Recently, calculation of 
Beenakker et. al. [p2| , p3| have solved the ID, N-channel case with amplification/absorption, 
both analytically and numerically, using the Dorokhov-Mello-Pereyra-Kumar (DMPK) |^ 
equation. Their study shows that for the N-channel case with absorption, reflection proba- 
bility distributions are Gaussian; and for the amplifying case reflection probability follow the 
Laguerre distribution. For = 1, these equations reduce to ours. All the theoretical studies 
25|^8| conclude that the amplification is enhanced by the localization effect of disorder. 



In conclusions, we have calculated the statistics of the reflection coefficient (r) and its 
associated phase (6) for a wave reflected from an amplifying/ absorbing disordered medium, 
for different disorder strengths and lengths of the sample. This numerical calculation clarifies 
several doubts recently raised on our previous paper on random phase approximation. 
We have modeled coherent amplification/absorption by adding a constant imaginary part 
to the real potential. We have discussed the cases of modeling with a complex potential. 
Following the Buttiker S-matrix approach, we have derived a Langevin equation for the 
reflection coefficient, which models absorption without the reflection and represent a case 
of pure stochastic absorption. Within the RPA, the FP equation for the case of stochastic 
absorption is the same as that for absorption modeled by a complex potential. Drawback of 
this numerical model calculation is that we can not go very large "r" due to the computation 
limitation for amplifying case, and for this reason we have not calculated the average "r" 

Our conclusions are the following: 
Anderson localization enhances the amplification /attenuation in a coherently amplify- 
ing/absorbing disordered media. The effective length of the pseudo cavity of the lasing action 
is the order of the localization length. 

Modeling active medium with a complex potential gives the following details results for 
both the amplifying and the absorbing cases: 

(a) Weak disorder and Weak active strength: 
Random-phase approximation is good for weak disorder {2k^ » 1). The strength of the 
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active part does not affect the random phase approximation up to a high value of the active 
parameter (D). Probabihty distribution P(r) of the reflection coefficient r has steady state 
solutions. 

(b) Weak disorder and Strong active strength: 

The random phase approximation (RPA) still holds for this region {2kp^ >> 1). Strength of 
the active part does not significantly affect the random phase approximation. P(r) saturates 
very fast. Wave gets absorbed before the reflection. 

(c) Strong disorder and Weak active strength: 

Random phase approximation does not hold in this regime. The strength of the active 
parameter does not affect much the phase distribution, which remains similar to that for the 
passive medium. There is a competition between amplification/absorption and localization 
which is shown by the double peaked distributions of P(r). 

(d) Strong disorder and Strong active strength: 

The probability distribution of the reflection coefficient peaks near r=l, for large lengths. 
The steady state reflection distribution {P{r)) is basically a delta function at r = 1 and the 
steady state phase distribution P{9) is a delta function at = vr, i.e., total reflection with 
opposite phase. The medium behaves as a perfect reflector. 

(e) Stochastic absorption: 

Modeling absorption without reflection shows that for the weak disorder case, the model is 
same as that of modeling absorption by adding an imaginary part to the real potential. 

(i) In case of strong disorder, and weak absorption, phase distribution P{9) does not 
change relative to the passive disordered medium and the probability distribution of the re- 
flection coefficient is similar to the model with the complex potential. The phase distribution 
typically has double peaks similar to that for a passive medium [Q . 

(ii) For strong disorder and strong absorption, the wave cannot penetrate deep inside 
the sample; it gets absorbed before it is reflected back. P(r) peaks near r = and the 
phase-distribution tries to peak far from 6 = n, which is different from modeling absorption 
by introducing a complex potential. 
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Finally, it would be interesting to study the problem in higher dimension with nonlinear 
amplification, we are looking for this. This problem also demand more experiments to look 
for. As fine optical fibers are easily available, it may not be difficult to do the experiment 
at lower dimensions. 



APPENDIX A: STOCHASTIC ABSORPTION 



The main idea is to simulate absorption by enlarging the S-matrix to include some "fake" , 
or "side" channels that remove few probability flux out of reckoning. Consider scattering 
channels 1 and 2 connected through current leads to two quantum "fake" channels 3 and 
4 that carry electrons to the reservoir, ( with chemical potential fi) as shown in the Fig.|l4 



(for only one scatterer). This is a phenomenological way of modeling absorption |]39| , ^ 
Electrons entering into the channels 3 and 4 are absorbed regardless of their phase and 
energy. The absorption is proportional to the strength of the "coupling parameter" e. A 
symmetric scattering matrix 5* for such a system can be written as : 



(Al) 



aR 


aT 







aT 


aR 


/5 










-aR* 


-aT* 


/5 





-aT* 


-aR* 



where R and T are the reflection and the transmissions amplitudes for the single scatterer 



, a = y/1 — e is the absorption coefficient and (3 = \fe. 

Now, the full scattering matrix is unitary for all positive real e < 1. 



But the sub-matrix connecting directly the channels 1 and 2 only is not unitary. We 
will now explore this fact. We observe from the above sub- matrix that for every scattering 
involving channels 1 and 2 only the reflection and the transmission amplitudes get multiplied 
by a (the absorption parameter). Now keeping this fact in mind, we will derive the Langevin 
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equation for the reflection amplitude R{L) for N scatterer (of length L) each characterized 
by a random S-matrix of this type, with the same a value. 

Consider N random scatterers for a ID sample of length L with (= Na, a = unitspacing) 
the reflection amphtude R{L), and A^" + 1 scatterers of the sample length L + AL, with the 
reflection amphtude R{L + AL). That is, let us start with the N scatterers, and add one 
more scatterer to the right to make up the N+1 scatterers. Now, we want to see the relation 
between R{L) and R{L + AL). That is, let there be a delta- function potential scatterer 
between L and L + AL, positioned at the point L + AL/2, that can be considered as the 
efi^ective scatterer due to the extra added scatterers for length AL. For kAL -C 1, we can 
treat the extra added scatterer as an effective delta potential Vo{L)S{x — L — AL/2) with 
Vo{L) — V{L)AL. (We consider the continuum hmit, a — > 0, A?" — > oo, Na — L, flxed ). 

Now, for a plane- wave scattering problem for a delta function potential of strength vq, 
which is at x = and has complex reflection and transmission amplitudes ro and to re- 
spectively, we have from the continuity condition for the wave function and discontinuity 
condition for the derivative of the wave function (which one gets by integrating Schrodinger 
equation across the delta function): 

and to = Tq + 1. 

Considering vq — v{L)AL the smallness parameter, one gets expressions up to flrst order in 
AL for To and to'. 

ro = ^ and 

^0 = 1 + where we have taken /2m — 1. 

Now to introduce absorption we will write: 
r ^ ret = ro(l - e)^/^, 
t^ta^ to{l - e)^/^ 

This means that for every scattering the reflection and transmission amplitudes get 



23 



modulated by a factor of a. 

Now consider a plane wave incident on the right side of the sample of length L + AL. 
Summing all the processes of direct and multiple reflections and transmissions, on the right 
side of the sample of length L, with the effective delta potential at L + AL/2, one gets, 

R[L + AL) = re^^^ 

+ (A2) 

Summing the above geometric series, substituting the values of r and t, and taking the 
continuum limit for L , one gets from the Langevin equation: 



= -aR{L) + 2ikR{L) + z|(r^(L))(l + R{L)f , (A3) 

with the initial condition R{L) = for L = 0, and a is the absorbing parameter. This is 
not quite the same as the Eq(4.4) obtained by introducing an imaginary potential. It turns 
out, however, this Langevin equation gives same results in the regime of weak disorder but 
differs qualitatively in the regime of strong disorder. 
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FIGURES 




FIG. 1. Coherently amplifying(marked by shade) medium in presence of disorder. Here we 
emphasize that the coherent-back-scattering (marked by the loop), hence localization will not be 
affected due to the presence of coherent amplification. As a result reflected wave will be amplified. 



28 




a b 
Wavelength — ^ X I2n (= 1/k) 

FIG. 2. A schematic picture of the localization window for a 3D Maxwell wave in {k, /e)-space, 
where A; = 27r/A is the incoming wave vector and Zg is the free scattering mean free path . Curve 
1 is for weak disorder, never satisfy the Mott-Iofee-Regal (MIR) criterion, i.e kle ^ 1. Curve 2 is 
for strong disorder case and has a localization window. The two mobility edges are denoted by a 
and h. The point a on the A axis is the correlation length of the disorder. For \ « a (geometric 
optics limit) the MIR condition for the localization, i.e., kle ^ Ij is not satisfied. For A > A;, (the 
Rayleigh limit) again the MIR criterion is not satisfied due to the weak Rayleigh scattering. 
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Non-radiative 
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FIG. 3. Proposal for a lasing scheme using Anderson localization for a non-resonant feedback 
mechanism: The figure shows a typical 3-level lasing scheme. Proposals for the lasing action are 
the following: Keep the pumping wave length Ap, which will excite electrons from the ground state 
to the excited state, inside the geometrical optic regime. The wave length of the lasing signal A^;, 
should lie inside the localization window, as shown in Fig.|2| 
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Partially 

Reflecting mirror transmitting mirror 

/ \ 



Amplifying medium 



Resonant cavity 



FIG. 4. Schematic picture of a Fabry-Perot etalon, a resonant cavity for the usual resonant 
feedback of lasing media. The ampHfying medium is kept between the two parallel mirrors which 
are used as resonators. 
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(a) 



Repeated Stimulated Emission 




A, ' ► A, 



Initial intensity distribution Gain distribution 

/V A- 4 

KX) l(A.)g(A.) l(A.)g"(A.) 

Gain narrowing through sucessive encounters 

FIG. 5. (a) A schematic picture in case of non-resonant "feedback" by feed-forward in a ho- 
mogeneous ampUfying medium. A wave is coherently amphfied as it moves forward. The dotted 
arrows denote the starting wave amplitude and the solid arrows show the coherent addition to 
the wave amplitude by repeated stimulated emission, (b) A typical gain narrowing process by 
successive coherent amplifications. 
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FIG. 6. Limiting probability distribution Pi^^ of the reflection coefHcient r, within the 

random phase approximation, for the cases (a) Coherent absorption (D = 1), (b) Unitary (D = 0), 
and (c) Coherent amplification {D = —1). 
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Phase 

FIG. 7. Probability distribution P{9) against the active parameter strength D in the weak 

disorder regime for a fixed disorder parameter strength 2kFS, = 50 and fixed sample length 1 = 1. 
Plots are for (a) Coherent amplification (D are -ve), and (b) Coherent absorption (D are +ve). 
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FIG. 8. Probability distribution P{r) against the sample length I, with fixed weak disorder 

strength 2kF(^ = 100 and fixed weak active parameter \D\. Plots are for (a) Coherent amplification 
(D = -0.5), and (b) Coherent absorption (D = +1). 
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FIG. 9. Probability distribution P(r) against the sample length Z, with fixed weak disorder 

strength 2kF^ = 100 and fixed strong active parameter \D\ = 10. Plots are for the cases (a) 



Coherent amplification (D = -10), and (b) Coherent absorption (D = +10). 
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FIG. 10. Probability distribution P(r) and P{9) separately against the sample length I. Pa- 
rameters are fixed strong disorder strength 2kpS^ = .001 and fixed weak active parameter \D\ = 0.5. 
Plots are for the cases (a) Coherent amplification {D = —0.5), and (b) Coherent absorption 
(D = +0.5). 
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FIG. 11. Probability distribution P{r) and P{6) separately against the sample length I. 

Parameters are fixed strong disorder strength = .001 and fixed strong active parameter 

\D\ = 10. Plots are for the cases (a) Coherent amplification {D = —10), and (b) Coherent 
absorption {D = +10) 
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FIG. 12. Probability distributions (a) P(r) and (b) P{6) against the sample length I with fixed 
strong disorder strength = .001 and fixed weak stochastic absorption parameter D = 0.5. 
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FIG. 13. Probability distributions (a) P(r) and (b) P{9) against the sample length I with fixed 
strong disorder strength = .001 and fixed strong stochastic absorption parameter D = 10. 
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FIG. 14. Modeling "absorption" by fake channels: Channels 1 and 2 are coupled through the 
current leads to two "fake" channels 3 and 4 which connect to a thermal reservoir. 
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